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O ■ 1 Introduction 

(N 

In spite of what the title might suggest, I shall not try to cover in these lec- 
tures all interesting aspects of the theory of the quark-gluon plasma. I shall 
rather focus on progress made in recent years in understanding the high tem- 
perature phase of QCD by using weak coupling techniques. Such techniques 
go far beyond strict perturbation theory viewed as an expansion in powers of 
the gauge coupling. In fact such an expansion becomes meaningless as soon 
^ i as the coupling is not vanishingly small. However, we shall see that a rather 

simple structure emerges from weak coupling studies, with a characteristic 
hierarchy of scales and degrees of freedom. The interactions renormalize the 
[~^ ' properties of these elementary degrees of freedom, but does not destroy the 

^D ■ simple picture of the high temperature quark-gluon plasma as a system of 

weakly interacting quasiparticles. As we shall see at the end of these lectures, 

this picture is supported by a first principle calculation of the entropy which 

T^ ' reproduces accurately lattice data above 2 or 3 times the critical temperature. 

'pH. Some of the material presented here is borrowed from the recent review ||l| , 

^ i' and complements can also be found in y,§,§J5[|2l. Another perspective on 

O . some of the topics discussed here can be found in the lectures by A. Rebhan. 

The outline of the lectures is the following. In order to get a first rough 
picture of the phase diagram of hadronic matter I use the bag model to 
describe the quark-hadron phase transition: this exercise will give us some 
familiarity with the thermodynamics of massless, non-interacting, particles. 
C^ ' Then I briefly recall some techniques of quantum field theory at finite tem- 

perature needed to treat the interactions [[7[|q,|9|, |lO|jll|Jl^ , and introduce the 
concept of effective theory in a simple case of a scalar field. Then I proceed 
to an analysis of the various important scales and degrees of freedom of the 
quark-gluon plasma and focus on the effective theory for the collective modes 
which develop at the particular momentum scale gT, where g is the gauge 
coupling and T the temperature. A powerful technique to construct the ef- 
fective theory is based on kinetic equations which govern the dynamics of the 
hard degrees of freedom. Some of the collective phenomena that are described 
by this effective theory are briefly mentioned. Then I turn to the calculation 
of the entropy and show how the information coded in the effective theory 
can be exploited in (approximately) self-consistent calculations p^,p^,p^ . 
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2 The quark-hadron transition in the bag model. 

The phase diagram of dense hadronic matter has the expected shape indi- 
cated in Fig. |l|. There is a low density, low temperature region, corresponding 
to the world of ordinary hadrons, and a high density, high temperature region, 
where the dominant degrees of freedom are quarks and gluons. The precise 
determination of the transition line requires elaborate non perturbative tech- 
niques, such as those of lattice gauge theories (see the lectures by F. Karsch). 
But one can get rough orders of magnitude for the transition temperature and 
density using a simple model dealing mostly with non-interacting particles 
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Fig. 1. The expected phase diagram of hot and dense hadronic matter in the plane 
{fj,B,T), where T is the temperature and /is the baryon chemical potential 



Let us first consider the transition in the case where ^b = 0. At low 
temperature this baryon free matter is composed of the lightest mesons, i.e. 
mostly the pious. At sufficiently high temperature one should also take into 
account heavier mesons, but in the present discussion this is an inessential 
complication. We shall even make a further approximation by treating the 
pion as a massless particle. At very high temperature, we shall consider that 
hadronic matter is composed only of quarks and antiquarks (in equal num- 
bers), and gluons, forming a quark-gluon plasma. In both the high tempera- 
ture and the low temperature phases, interactions are neglected (except for 
the bag constant to be introduced below). The description of the transition 
will therefore be dominated by entropy considerations, i.e. by counting the 
degrees of freedom. 

The energy density e and the pressure P of a gas of massless pions are 
given by: 

30 ' 90 ' ^ ' 

where the factors 3 account for the 3 types of pions (tt"*", tt", w^). 

The energy density and pressure of the quark-gluon plasma are given by 
similar formulae: 
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S7.-T^-B, 



(2) 



where 37 = 2x8 + |x2x2x2x3is the effective number of degrees of freedom 
of gluons (8 colors, 2 spin states) and quarks (3 colors, 2 spins, 2 flavors, q 
and q). The quantity B, which is added to the energy density, and subtracted 
from the pressure, summarizes interaction effects which are responsible for a 
change in the vacuum structure between the low temperature and the high 
temperature phases. It was introduced first in the "bag model" of hadron 
structure as a restoring force needed to equilibrate the pressure generated by 
the kinetic energy of the quarks inside the bag IM]. Roughly, the energy of 
the bag is 



E{R)^fR^B + ^, 



(3) 



where C/R is the kinetic energy of massless quarks. Minimizing with respect 
to R, one finds that the energy at equilibrium is E (Rq) = 4i?Vo, where 
Vq — 47ri?Q/3 is the equilibrium volume. For a proton with i?o ~ 1 GeV and 
Rq « 0.7 fm, one finds Eq/Vq ~ 0.7 GeV/fm^, which corresponds to a "bag 
constant" B « 175 MeV/fm^, or B^/"^ « 192 MeV. 

We can now compare the two phases as a function of the temperature. 
Fig. H shows how P varies as a function of T^. One sees that there exists a 
transition temperature 



Tr 



\ 17^2 ) 



1/4 



Si/4 w 0.72 B^/\ 



(4) 



beyond which the quark-gluon plasma is thermodynamically favored (has 
largest pressure) compared to the pion gas. For B^'^ w 200 MeV, Tc ~ 150 
MeV. 



quark-gluon 
plasma 




Fig. 2. The pressure of the massless pion gas compared to that of a quark-gluon 
plasma, showing the transition temperature T^. 
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The variation of the entropy density s = dP/dT as a function of the tem- 
perature is displayed in Fig. 0. Note that the bag constant B does not enter 
expHcitly the expression of the entropy. However, B is involved in Fig. S indi- 
rectly, via the temperature Tc where the discontinuity As occurs. One verifies 
easily that the jump in entropy density As = Ae/Tc is directly proportional 
to the change in the number of active degrees of freedom when T crosses Tc. 

In order to extend these considerations to the case where ^b 7^ 0, we 
note that the transition is taking place when the total pressure approxi- 
mately vanishes, that is when the kinetic pressure of quarks and gluons ap- 
proximately equilibrates the bag pressure. Taking as a criterion for the phase 
transition the condition P = 0, one replaces the value (^ for Tc by the value 
(90/377r2)i/45i/4 p. OJOB^/'^, which is nearly identical to (|). We shall then 
assume that for any value of /is and T, the phase transition occurs when 
P{lJ,B,T) — B, where B is the bag constant and P{^b,T) is the kinetic 
pressure of quarks and gluons: 

P(MB,T).|.^T^ + f(T^-,iJ). (5) 

The transition line is then given by P{p-c, Tc) = B, and it has indeed the 
shape illustrated in Fig. [|. 

The model that we have just described reproduces some of the bulk fea- 
tures of the equation of state obtained through lattice gauge calculations (see 
the lecture by F. Karsch). In particular, it exhibits the characteristic increase 
of the entropy density at the transition which corresponds to the emergence 
of a large number of new degrees of freedom associated with quarks and glu- 
ons. Its simplicity has made it popular for instance among the practitioners 
of hydrodynamic calculations with which one tries to simulate the behavior of 
matter produced in high energy nuclear collisions. As such it has been very 
useful. One should be cautious however when attempting to draw too de- 
tailed conclusions about the nature of the phase transitions from such simple 
models. In particular this model predicts (by construction!) a discontinuous 
transition; but this prediction should not be trusted. Further discussion of 
this model can be found in B 



3 Quantum Fields at Finite Temperature 

The effects of interactions among quarks and gluons at finite temperature can 
be calculated by using the tools of quantum field theory at finite temperature. 
We shall briefly recall some essential formalism, and emphasize in particular 
the periodicity properties of the propagators. At the end of this section we 
discuss, with a simple example of a scalar field, the method of effective field 
theory which proves useful in problems where various scales can be separated. 
In the example that we shall consider, the separation of scale is provided by 
the Matsubara frequencies. As we shall see, in some cases, one is lead to 
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Fig. 3. The entropy density. The jump As at the transition is proportional to the 
increase in the number of active degrees of freedom 



single out the mode with vanishing Matsubara frequency. The corresponding 
effective theory is a classical field in three dimensions, and the procedure 
commonly called 'dimensional reduction'. 



3.1 Finite Temperature calculations 

All thermodynamic observables can be deduced from the partition function: 
Z = tre~'^^. (6) 

Thus the energy density and the pressure are given by: 






1 d 



(7) 



In order to calculate the partition function, one may observe that e ^^ is 
like an evolution operator in imaginary time: 



t 



'i(i 



-iHt 



-f3H 



(8) 



One may then take advantage of all the techniques developed to evaluate 
matrix elements of the evolution operator in quantum mechanics or field 
theory. 

For instance one may use a perturbative expansion. We assume that one 
can split the hamiltonian into H — H^ + Hi with Hi <C Hq , and define the 
following "interaction representation" of the perturbation Hi : 



Hi{r) 



e'^^^Hie-^"" 



(9) 



and similarly for other operators. Using standard techniques, one can then 
obtain the following expression for the partition function Z: 



Zo(TexpJ- / dTiJi(T) l)o. 



(10) 
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In this equation, the symbol T imphes an ordering of the operators on its 
right, from left to right in decreasing order of their time arguments; Zq = 
tr e~^^° and, for any operator O, 



{O)o ^ Tr [^^Oj . (11) 

One commonly refers to r as the "imaginary time" (r is real). This r has 
no direct physical interpretation: its role here is to properly keep track of 
ordering of operators in the perturbative expansion. 

In field theory, it is often more convenient to use the formalism of path 
integrals. Let us recall for instance that for one particle in one dimension the 
matrix element of the evolution operator can be written as 

{q2\e-'"'\qi) = r ''v{q{t)) e*L?(^"9^-n?))dt ^ (12) 

where qi and <?2 denote the positions of the particle at times and t respec- 
tively. Changing t — > — ir, and taking the trace, one obtains the following 
formula for the partition function: 

Z = tr e-P" = / V{q) exp J - / f ^mg^ + v{q)] \ . (13) 

Jq{l3)=q{0) [ Jo \^ J) 

This expression immediately generalizes to the case of a scalar field, for 
which the Lagrangian is of the form: 

1 777'^ 



1,„ .., 1 



-{d,<i>f - -{Vct>f - —4>'' -V{^). (14) 

Again, we replace t by —it, Oq — dt by idr^ so that {do<j))^ -^ —{dr(j))^- The 
partition function becomes then (integrations over spatial coordinates are 
implicit): 

(15) 
where the integral is over periodic fields: 0(0) = 4'{0)- 



Remarks, i) The partition function (|15|) may be viewed formally as a sum 

over classical field configurations in four dimensions, with particular boundary 

conditions in the (imaginary) time direction. 

ii) At high temperature, /3 — > 0, the time dependence of the fields play no 

role. The partition function becomes that of a classical field theory in three 

dimensions: 

2 = |o(0)exp|-/3|dVQ(V<^)V^0^ + V^(0)j|. (16) 



Theory of the quark-gluon plasma 7 

Ignoring the time dependence of the fields amounts to take into account only 
the Matsubara frequency iuj,^ = 0. We shall discuss later explicit examples of 
this "dimensional reduction". 

iii) Note the Euclidean metric in ( |l5[ ) . Since the integrand is the exponential of 
a negative definite quantity, it is well suited to numerical evaluations, using for 
instance the lattice technique. 



3.2 Free propagators 

An important feature of the path integral representation of the partition 
function is the boundary conditions to be imposed on the fields over which 
one integrates. For the scalar case considered here, the field has to be periodic 
in imaginary time, with a period /3. Similar conditions hold for the fermion 
fields, which are antiperiodic in imaginary time, with the same period /3. It 
is instructive to see how these periodicity conditions emerge in the operator 
formalism, and for this reason we consider now the free propagators, first in 
the simple case of the non relativistic many body problem. The generalization 
to relativistic field is straightforward. 

Let us consider a system with unperturbed hamiltonian: 

Ho = '^ekalak, (17) 

fe 

where k denotes the set of quantum numbers necessary to specify a single 
particle state, for instance the three components of the momentum. We define 
time dependent creation and annihilation operators in the interaction picture: 



akir) = e"^«afce-^^° = e-'-^ak- (18) 

The last equalities follow (for example) from the commutation relations: 

[Ho,al] = Ekal [Ho,ak] == -efca^ (19) 

which hold for bosons and fermions. The single particle propagator can then 
be obtained by a direct calculation: 

Gfe(Ti - T2) = (Tafc(ri)aj,(T2))o 

= e-'^^^'-^^^ [0{ti - T2)(l ± rife) ± nkeiT2 - n)] , (20) 

where: 

Uk = (4afe)o = ^^,^ ^ ^ , (21) 

and the upper (lower) sign is for bosons (fermions). One can verify on the 
expression ( po| ) that, in the interval — /? < t = ri — T2 < /3, Gu(j) is a periodic 
(boson) or antiperiodic (fermion) function of r: 

Gfc(T-/3) = ±Gfe(r) (0<T</3). (22) 
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(To show this relation note that e^'^^Uk = 1 ± n^.) It can therefore be repre- 
sented by a Fourier series 

Gk{r)^'^Y.^-'^"^Gk{ioj,), (23) 

where the uj^^s are cahed the Matsubara frequencies: 



uji, ~ 2vn / (3 bosons, 

Wjy — {2v + 1)tt/P fermions. 

The inverse transform is given by 

/3 



(24) 



G{iLu,)= dTe'--^G(T) = - — . (25) 

Using the property 

'^M^^E^"""' -/3<T</3 (26) 

1/ 

and (p3|), it is easily seen that G{t) satisfies the differential equation 

{dr + Ho)G{T)^6{T), (27) 



which may be also verified directly from (EQ). Alternatively, the single prop- 
agator at finite temperature may be obtained as the solution of this equation 
with periodic (bosons) or antiperiodic (fermions) boundary conditions. 

Remark. The periodicity or antiperiodicity that we have uncovered on the 
explicit form of the unperturbed propagator is, in fact, a general property of the 
propagators of a many-body system in thermal equilibrium. It is a consequence 
of the commutation relations of the creation and annihilation operators and the 
cyclic invariance of the trace. 

The propagator of the free scalar field A(t) = (T0(ri)(/i(r2), where t = 
Ti — T2 satisfies the differential equation 

[-dl - Vl+m^] Z\(riri;T2r2) = <5(ti - r2)^(ri - ra), (28) 

and obeys periodic boundary conditions. It admits the Fourier representation 

n 

where tUn — 27rn//3 and 

AiiLu„) = , ^ , . (30) 

By inverting the Fourier transform (pO[), one gets 

A{t) = ^ |(1 + iVfc)e-^''l^l -I- iVfee^"!^!} , (31) 

with Nk = l/(e'3^'' - 1). 
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3.3 Classical field approximation and dimensional reduction 

In the high temperature limit, /? — > 0, the imaginary-time dependence of the 
fields frequently becomes unimportant and can be ignored in a first approx- 
imation. The integration over imaginary time becomes then trivial and the 
partition function (|l3) reduces to: 



Z 



M fv{(t)) exp I -/3 / d^xHix) \ , (32) 



where 4> = </>(x) is now a three-dimensional field, and 



1 777 

n = -{V^r + — cj^' + Vicf>). (33) 

The functional integral in ( p2| ) is recognized as the partition function for 
static three-dimensional field configurations with energy / d^xHix). We shall 
refer to this limit as the classical field approximation. 

Ignoring the time dependence of the fields is equivalent to retaining only 
the zero Matsubara frequency in their Fourier decomposition. Then the Fourier 
transform of the free propagator is simply: 

Go(k) = ^ . (34) 



This may be obtained directly from (|29|) keeping only the term with a;^ = 0, 
or from eq. ( |3l| ) by ignoring the time dependence and using the approximation 

iV(s.) = ^-J. (35) 

Both approximations make sense only for e^ <^ T, implying N{ek) 3> 1. In 
this limit, the energy per mode is oc ekN{ek) ~ T, as expected from the 
classical equipartition theorem. 

The classical field approximation may be viewed as the leading term in 
a systematic expansion. To see that, let us expand the field variables in the 
path integral dlq) in terms of their Fourier components: 

^ir)=^^J2^-''^"'^(iiu,), (36) 

where the uj^^s are the Matsubara frequencies. The path integral ( p^ ) can 
then be written as: 



Z = M / -Dicbo) exp {-S[cbo]} , (37) 

where 4>o = 4>{uJi, = 0) depends only on spatial coordinates, and 

'xCe{x)\ . (38) 



exp{-S[(f>o]} = N2 /p(0^^o)expi - / At h?x 
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The quantity S[(j)o] may be called the effective action for the "zero mode" 
^0- Aside from the direct classical field contribution that we have already 
considered, this effective action receives also contributions from integrating 
out the non- vanishing Matsubara frequencies. Diagrammatically, ^[^o] is the 
sum of all the connected diagrams with external lines associated to 0Oi and 
in which the internal lines are the propagators of the non-static modes (/iiy^o ■ 
Thus, a priori, S[4>o] contains operators of arbitrarily high order in 0o, which 
are also non-local. In practice, however, one wishes to expand S[(t)o\ in terms 
of local operators, i.e., operators with the schematic structure am^„V™0o 
with coefficients am, n to be computed in perturbation theory. 

To implement this strategy, it is useful to introduce an intermediate scale 
A (A <^ T) which separates hard [k ^ A) and soft {k ^ A) momenta. All 
the non-static modes, as well as the static ones with k ^ A are hard (since 
K^ = ojI + k^ ^ A^ for these modes), while the static {uj^ = 0) modes with 
k ^ A are soft. Thus, strictly speaking, in the construction of the effective 
theory along the lines indicated above, one has to integrate out also the 
static modes with k ^ A. The benefits of this separation of scales are that 
(a) the resulting effective action for the soft fields can be made local (since the 
initially non-local amplitudes can be expanded out in powers of p/K, where 
p <C ^ is a typical external momentum, and K <^ A is a hard momentum 
on an internal line), and (6) the effective theory is now used exclusively at 
soft momenta, where classical approximations such as (|3^) are expected to be 
valid. This strategy, which consists in integrating out the non-static modes in 
perturbation theory in order to obtain an effective three-dimensional theory 
for the soft static modes (with cu,^ = Q and fc = |k| ^ yl), is generally referred 
to as "dimensional reduction" 0,|l|,^9|j2^,|l|,|2) . 

As an illustration let us consider a massless scalar theory with quartic 
interactions; that is, to = and V{(l)) — (.g^/4 !)(/>'* in (p^). The ensuing 
effective action for the soft fields (which we shall still denote as 0o) reads 

5[0o] = P:f{A) 

(39) 

where J-'{A) is the contribution of the hard modes to the free-energy, and AC 
contains all the other local operators which are invariant under rotations and 
under the symmetry —^ —0, i.e., all the local operators which are consistent 
with the symmetries of the original Lagrangian. We have changed the normal- 
ization of the field {(pa — > VT4>o) with respect to (p^)~(p3|), so as to absorb 
the factor (3 in front of the effective action. The effective "coupling constants" 



in (39), i.e. M^(yl), g'^{A), h{A) and the infinitely many parameters in AC, 
are computed in perturbation theory, and depend upon the separation scale 
A, the temperature T and the original coupling g^. To lowest order in g, 
gl « g^T, ft, « (the first contribution to h arises at order g^, via one- loop 
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diagrams), and M ^ gT, as we shall see shortly. Note that eq. ( p9| ) involves in 
general non-renormaUzable operators, via AC This is not a difficulty, how- 
ever, since this is only an effective theory, in which the scale A acts as an 
explicit ultraviolet (UV) cutoff for the loop integrals. Since however the scale 
A is arbitrary, the dependence on A coming from such soft loops must cancel 
against the dependence on A of the parameters in the effective action. 




One-loop tadpole diagram for the self-energy of the scalar field 



Let us verify this cancellation explicitly in the case of the thermal mass 
M of the scalar field, and to lowest order in perturbation theory. To this 
order, the scalar self-energy is given by the tadpole diagram in Fig. H. The 
mass parameter M^ (A) in the effective action is obtained by integrating over 
hard momenta within the loop in Fig. B[ 

M^A) 



2 
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- A)S^o 


2 

9' 
2 


(2^)3 

(Nik) 


u 

1 

+ 2fc 
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- e{A- 


T 



(40) 

where the 0-function in the second line has been generated by writing 0{k — 
A) = 1 ~ 0{A — k). The first term, involving the thermal distribution, gives 
the contribution 



M^ ^ '^ ^L^tl^ ^ y_T\ (41) 



g2 f d^fc N{k) _ g^ 
(27r)3 k ~ 24 

As it will turn out, this is the leading-order (LO) scalar thermal mass, and 
also the simplest example of what will be called "hard thermal loops" (HTL). 
The second term, involving 1/2/c, in ( ^0| ) is quadratically UV divergent, but 
independent of the temperature; the standard renormalization procedure at 
T = amounts to simply removing this term. The third term, involving the 
^-function, is easily evaluated. One finally gets: 

47r2 ~ 24 



M\A) ^M' -^AT^^-—{l---\. (42) 



The yl-dependent term above is subleading, by a factor A/T ^ 1. 

The one-loop correction to the thermal mass within the effective theory is 
given by the same diagram in Fig. H, but where the internal field is static and 
soft, with the massive propagator 1/(A:^ -I- M^(yl)), and couphng constant 
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33 ~ 9^T. Since the typical momenta in the integral will be fc ^ M, and 
M ^ M r^ gT, we choose A > gT. We then obtain 



2 '^ d3/fc „, . ., T 



SM^A) = ^ / j^ 0{A - k) 



2 J (27r)3 ' ' B+A'P{A) 

g^TA f nM A \ g^TA g^ 



, ,1 r arctan — ~ - — ^^MT , (43) 

47r2 \ 2A Mj 47r2 Stt ' ^ ' 

where the terms neglected in the last step are of higher order in AI/A or 
A/T. 

As anticipated, the yl-dependent terms cancel in the sum Af^ = M'^{A) + 
6M^{A), which then provides the physical thermal mass within the present 
accuracy: 

M^ = M^(A) + SM'^(A) ^ '^ ^MT. (44) 

^ '^ ^ ^ 24 Stt ^ ^ 

The LO term, of order g'^T'^, is the HTL M. The next-to-leading order (NLO) 
term, which involves the resummation of the thermal mass M(A) in the soft 
propagator, is of order g^MT ~ g^T^, and therefore non-analytic in g^. This 
non-analyticity is related to the fact that the integrand in (|4|) cannot be 
expanded in powers of M"^ /k'^ without generating infrared divergences. 

4 Effective theories for the quark-gluon plasma 

We return now to the quark-gluon plasma and analyze the various scales 
and degrees of freedom which are relevant in the weak coupling regime. We 
show that there is a hierarchy of scales controlled by powers of the gauge 
coupling g. We focus in these lectures on two particular momentum scales, 
the 'hard' one which is that of the plasma particles with momenta fc ~ T, 
and the 'soft' one with k ~ gT at which collective phenomena develop. We 
shall be in particular interested in the effective theory obtained when the hard 
degrees of freedom are 'integrated out'. The resulting effective theory describe 
long wavelength, low frequency collective phenomena; that is, it accounts for 
time dependent fields, in contrast to the example discussed in the previous 
section which concerned only static fields. As we shall see later, getting a 
complete description of the dynamics of the collective excitations turns out 
to be important also for the calculation of the equilibrium properties of the 
quark-gluon plasma. 

4.1 Scales and degrees of freedom in ultrarelativistic plasmas 

A property of QCD which is essential in the present discussion is that of 
asymptotic freedom, according to which the coupling constant depends on 
the scale p. as 

_ _ c,2 ^ 

47r \\i[fi/AQCD)' 



Theory of the quark-gluon plasma 13 

At high temperature, the natural scale is p, = 2t:T, so that the coupling 
becomes weak when 2'kT 3> Aqcd- At extremely high temperature the in- 
teractions become negligible and hadronic matter turns into an ideal gas of 
quarks and gluons: this is the quark-gluon plasma. As we shall see an impor- 
tant effect of the interactions is to turn free quarks and gluons into weakly 
interacting quasiparticles. 

In the absence of interactions, the plasma particles are distributed in mo- 
mentum space according to the Bose-Einstein or Fermi-Dirac distributions: 

Nk = ^5; -, nk = -5— -— r, (46) 

where Sk — k = |k| (massless particles), (3 = 1/T, and chemical potentials 
are assumed to vanish. In such a system, the particle density n is determined 
by the temperature: n (x T^. Accordingly, the mean interparticle distance 
n~^^^ ^ 1/T is of the same order as the thermal wavelength At = 1/fc of 
a typical particle in the thermal bath for which k ^ T. Thus the particles 
of an ultrarelativistic plasma are quantum degrees of freedom for which in 
particular the Pauli principle can never be ignored. 

In the weak coupling regime {g <C 1), the interactions do not alter sig- 
nificantly the picture. The hard degrees of freedom, i.e. the plasma particles, 
remain the dominant degrees of freedom and since the coupling to gauge 
fields occurs typically through covariant derivatives, Dx = dx + igA{x), the 
effect of interactions on particle motion is a small perturbation unless the 
fields are very large, i.e., unless A r^ T/g, where g is the gauge coupling: 
only then do we have dx ^ T ^ gA, where dx is a space-time gradient. We 
should note here that we rely on considerations, based on the magnitude of 
the gauge fields, which depend on the choice of a gauge. What is meant is 
that there exists a large class of gauge choices for which they are valid. And 
we shall verify a posteriori that within such a class, the final results are gauge 
invariant. 

Considering now more generally the effects of the interactions, we note 
that these depend both on the strength of the gauge fields and on the wave- 
length of the modes under study. A measure of the strength of the gauge 
fields in typical situations is obtained from the magnitude of their thermal 
fiuctuations, that is A = ^y{A'^{t,K)). In equilibrium (A^(i,x)) is indepen- 
dent of t and X and given by {A^) — G{t = 0, x = 0) where G{t, x) is the 
gauge field propagator. In the non interacting case we have (with Sk — k): 

Here we shall use this formula also in the interacting case, assuming that the 
effects of the interactions can be accounted for simply by a change of Sk- We 
shall also ignore the (divergent) contribution of the vacuum fiuctuations (the 
term independent of the temperature in (p7|)). 
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For the plasma particles Ek — k ^ T and {A^)t ~ T^- The associated 
electric (or magnetic) field fluctuations are {E^)t ~ ((9^)^)t "^ k'^{A^)T ^ 
T'^ and are a dominant contribution to the plasma energy density. As already 
mentioned, these short wavelength, or hard, gauge field fluctuations produce 
a small perturbation on the motion of a plasma particle. However, this is 
not so for an excitation at the momentum scale k ~ gT, since then the two 
terms in the covariant derivative dx and gAx become comparable. That is, 
the properties of an excitation with momentum gT are expected to be non 
perturbatively renormalized by the hard thermal fluctuations. And indeed, 
the scale gT is that at which collective phenomena develop. The emergence 
of the Debye screening mass rriD ^ gT is one of the simplest examples of 
such phenomena. 

Let us now consider the fluctuations at this scale gT <C T, to be re- 
ferred to as the soft scale. These fluctuations can be accurately described by 
classical fields. In fact the associated occupation numbers Nk are large, and 
accordingly one can replace Nk by T/sk in ( |47| ) . Introducing an upper cut-off 
gT in the momentum integral, one then gets: 



gT 



T 



{A')gT-J d'k^^gT'. (48) 

Thus AgT ~ \/gT so that gAgT ~ g'^^^T is still of higher order than the 
kinetic term dx ~ gT. In that sense the soft modes with k ^ gT are still per- 
turbative, i.e. their self-interactions can be ignored in a first approximation. 
Note however that they generate contributions to physical observables which 
are not analytic in g^ , as shown by the example of the order g'^ contribution 
to the energy density of the plasma: 

e(3) ^ r' d^fc uo^i _l ^ ^ ojIi LUpi — ^ g^T\ (49) 







'^P'e^.,/T_i-^'^P/'^P' — 



where lUpi ~ gT is the typical frequency of a collective mode. 

Moving down to a lower momentum scale, one meets the contribution of 
the unscreened magnetic fluctuations which play a dominant role for k ^ g^T. 
At that scale, to be referred to as the ultrasoft scale, it becomes necessary to 
distinguish the electric and the magnetic sectors (which provide comparable 
contributions at the scale gT). The electric fluctuations are damped by the 
Debye screening mass (el — k'^ + rnjj w rrijy when k ~ g^T) and their 
contribution is negligible, of order g*T^. However, because of the absence of 
static screening in the magnetic sector, we have here e/c ~ fc and 



(A^T^T/ d^fc-^g^T^, (50) 



so that gAg2rp ~ g'^T is now of the same order as the ultrasoft derivative 
dx ~ g^T: the fluctuations are no longer perturbative. This is the origin of 
the breakdown of perturbation theory in high temperature QCD. 
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Fig. 5. Example of a multiloop diagram which is infrared divergent 

To appreciate the difficulty from another perspective, let us first observe 
that the dominant contribution to the fluctuations at scale g^T comes from 
the zero Matsubara frequency: 

iA%.r = TY^f^^k-^^^Tf,^k^^. (51) 

Thus the fluctuations that we are discussing are those of a three dimensional 
theory of static fields. Following Linde |^,^ consider then the higher order 
corrections to the pressure in hot Yang-Mills theory. Because of the strong 
static fluctuations most of the diagrams of perturbation theory are infrared 
(IR) divergent. By power counting, the strongest IR divergences arise from 



ladder diagrams, like the one depicted in Fig. 4.1, in which all the propaga- 
tors are static, and the loop integrations are three-dimensional. Such n-loop 
diagrams can be estimated as (^ is an IR cutoff): 

which is of the order g^T^ \n{T/^i) if n = 4 and of the order g^T'^ {g'^T/^i) 



if n > 4. (The various factors in (|52| ) arise, respectively, from the 2(n — 1) 



three-gluon vertices, the n loop integrations, and the 3(n — 1) propagators.) 
According to this equation, if /i ~ g^T , all the diagrams with n > 4 loops 
contribute to the same order, namely to 0{g^). In other words, the correction 
of 0{g^) to the pressure cannot be computed in perturbation theory. 

4.2 Effective theory at scale gT 

Having identified the main scales and degrees of freedom, our task will be 
to construct appropriate effective theories at the various scales, obtained by 
eliminating the degrees of freedom at higher scales. We shall consider here 
the effective theory at the scale gT obtained by eliminating the hard degrees 
of freedom with momenta k ^ T. 

The soft excitations at the scale gT can be described in terms of average 
fields |^,g6[. Such average fields develop for example when the system is 
exposed to an external perturbation, such as an external electromagnetic 
current. In QED, we can summarize the effective theory for the soft modes 
by the equations of motion: 

d^F^''=j:n,+jLt (53) 
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that is, Maxwell equations with a source term composed of the external per- 
turbation Jej-t, and an extra contribution j'j^^j^ which we shall refer to as the 
induced current. The induced current is generated by the collective motion 
of the charged particles, i.e. the hard degrees of freedom. It may be regarded 
itself as a functional of the average gauge fields and, once this functional is 
known, the equations above constitute a closed system of equations for the 
soft fields. 

The main problem is to calculate jmd- This is done by considering the 
dynamics of the hard particles in the background of the soft fields. For QED, 
the induced current can be obtained using linear response theory. To be 
more specific, consider as an example a system of charged particles on which 
is acting a perturbation of the form J dxjfj,{x)A'^{x), where jfj.{x) is the 
current operator and A^{x) some applied gauge potential. Linear response 
theory leads to the following relation for the induced current: 

jl^'' = Jd^yn^.{x-y)A''{y), 

n^A^ -y) = -i0{xo - yo){[j,^{x),jAy)])eq., (54) 

where the (retarded) response function 77^^ {x — y) is also referred to as the 
polarization operator. Note that in (^^, the expectation value is taken in 
the equilibrium state. Thus, within linear response, the task of calculating 
the basic ingredients of the effective theory for soft modes reduces to that of 
calculating appropriate equilibrium correlation functions. 

In fact we shall need the response function only in the weak coupling 
regime, and for particular kinematic conditions which allow for important 
simplifications. In leading order in weak coupling, the polarization tensor is 
given by the one-loop approximation. In the kinematic regime of interest, 
where the incoming momentum is soft while the loop momentum is hard, we 
can write II{cu,p) = g'^T'^f{cu/p,p/T) with / a dimensionless function, and 
in leading order in p/T ~ g, 77 is of the form g^T"^ f (uj / p) . This particular 
contribution of the one-loop polarization tensor is an example of what has 
been called a "hard thermal loop" [|^J2U2|,|^,|l|j32[|2|,|6) ; for photons in 



QED, this is the only one. It turns out that this hard thermal loop can 
be obtained from simple kinetic theory, and the corresponding calculation is 
done in the next subsection. 

In non Abelian theory, linear response is not sufficient: constraints due to 
gauge symmetry force us to take into account specific non linear effects and 
a more complicated formalism needs to be worked out. Still, simple kinetic 
equations can be obtained in this case also, but in contrast to QED, the 
resulting induced current is a non linear functional of the gauge fields. As a 
result, it generates an infinite number of "hard thermal loops" . 
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5 Kinetic equations for the plasma particles 

The hard degrees of freedom enter the equations of motion (B3[) for the 
soft collective excitations only through their average density or current, that 
is, through the induced current. This induced current can be calculated by 
studying the dynamics of the plasma particles in the background of soft exter- 
nal gauge fields. This is what we now turn to. In order to keep the discussion 
at an elementary level, we shall merely analyze the main steps involved in 
the derivation of the corresponding QCD equations in the simpler context of 
non relativistic electromagnetic plasmas. The QCD equations are presented 
at the end of this section. 

5.1 One-loop polarization tensor from kinetic theory 

As indicated above, in the kinematic regime considered, the dominant contri- 
bution to the one loop polarization tensor can be obtained using elementary 
kinetic theory, and we present now this calculation. We consider an electro- 
magnetic plasma and momentarily assume that we can describe its charged 
particles in terms of classical distribution functions fq{p,x) giving the den- 
sity of particles of charge q {q = ±e) and momentum p at the space-time 
point X = (i,r) pj. We consider then the case where collisions among the 
charged particles can be neglected and where the only relevant interactions 
are those of particles with average electric (E) and magnetic (B) fields. Then 
the distribution functions obey the following simple kinetic equation, known 
as the Vlasov equation |g^ '■ 

where v = de^/dp is the velocity of a particle with momentum p and energy 
£p (for massless particles v — p), and F = q(Ei + v A B) is the Lorentz force. 
The average fields E and B depend themselves on the distribution functions 
fq. Indeed, the induced current 

where v^ = (l,v), is the source term in the Maxwell equations (^Sl) for the 
mean fields. 

When the plasma is in equilibrium, the distribution functions, denoted as 
fq(p) = f°{ep), are isotropic in momentum space and independent of space- 
time coordinates; the induced current vanishes, and so do the average fields 
E and B. When the plasma is weakly perturbed, the distribution functions 
deviate slightly from their equilibrium values, and we can write: /g(p, a;) = 
/''(sp) + Sfq{p,x). In the linear approximation, Sf obeys 

{v-d,)Sfq{p,x) = -q^-E^, (57) 



jU^) -^ I 7K^^' (/+(P' ^) - /-(P: ^)) > (56) 
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where z; • 9^ = 9t + v • V. The magnetic field does not contribute because of 
the isotropy of the equilibrium distribution function. 
It is convenient here to set 

df" 
Sfg{p,x) = -qW{x,^)^, (58) 

thereby introducing a notation which will be useful later for the QCD case. 
Since 

/,(P, x) = f{sp) ~ qW{x, v) ^ ^ fiep - qWix, v)), (59) 

W{x,'v) may be viewed as a local distortion of the momentum distribution 
of the plasma particles. The equation for W is simply: 

{v■^^)W{x,^^) = vE(a;). (60) 

Contrary to (^) , the linearized equations ( |57| ) or ( |60| ) do not involve the 
derivative of / with respect to p, and they can be solved by the method of 
characteristics: v-dx is the time derivative of J/(p, x) along the characteristic 
defined by dx/dt — v. Assuming then that the perturbation is introduced 
adiabatically so that the fields and the fluctuations vanish as e***'' (77 ^ 0+) 
when Iq — > — cx3, we obtain the retarded solution: 



(61) 



W{x,v) = I di'v-E(x-v(i-i')>^'): 
J —00 

and the corresponding induced current: 

r rl^r, dfO 1"°° 
jUx)--^^'Jj^^^i;l drv.E(.-.r). (62) 



Since E = — VA° — dA/dt, the induced current is a linear functional of A^. 
At this point we assume explicitly that the particles are massless. In this 
case, V is a unit vector, and the angular integral over the direction of v 
factorizes in (|6^). Then, using (g4|) as definition for the polarization tensor 
7T'"'(a; — y), and the fact that the Fourier transform of J^, dr e^^''^ f {x — vt) 
is i f{Q)/{v ■ Q + iri), with Q^ ~ (^iQ) a-nd f{Q) the Fourier transform of 
f{x), one gets, after a simple calculation |Q : 

nf,^{uj,q)^ml)\-Sf,QS^o+^ -. ^ " , . \ , (63) 

|_ J in uj — v ■ q + ir] ) 

where the angular integral J df2 runs over all the orientations of v, and tti^j 
is the Debye screening mass: 



.2 _ ^ r App^'^ 

TT^ ,/n de^ 
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It turns out that ( |63| ) is the dominant contribution at high temperature to 
the one-loop polarization tensor in QED, provided one substitutes for f'^ the 
actual quantum equilibrium distribution function, that is, f'^{sp) = rip, with 
np given in (Hq). After insertion in (p#, this yields to|, — e^T^/3. 

In the next subsection, we shall address the question of how simple kinetic 
equations emerge in the description of systems of quantum particles, and 
under which conditions such systems can be described by seemingly classical 
distribution functions where both positions and momenta are simultaneously 
specified. 

We shall later find that the expression obtained for the polarization tensor 
using simple kinetic theory generalizes to the non Abelian case. This is so 
in particular because the kinematic regime remains that of the linear Vlasov 
equation, with straight line characteristics. 

5.2 Kinetic equations for quantum particles 

In order to discuss in a simple setting how kinetic equations emerge in the 
description of collective motions of quantum particles, we consider in this 
subsection a system of non relativistic fermions coupled to classical gauge 
fields. Since we are dealing with a system of independent particles in an 
external field, all the information on the quantum many-body state is encoded 
in the one-body density matrix |9|JlO(| : 



p{v,v',t) = {W\v',m{v,t)), (65) 

where ^ and tf'T are the annihilation and creation operators, and the average 
is over the initial equilibrium state. It is on this object that we shall later 
implement the relevant kinematic approximations. To this aim, we introduce 
the Wigner transform of p(r, r', t) pq,pq|: 

/(p, R, t) = I d^s e-'P-^ p (r + |, R - £, i) . (66) 

The Wigner function has many properties that one expects of a classical phase 
space distribution function as may be seen by calculating the expectation 
values of simple one-body observables. For instance the average density of 
particles n(R) is given by: 

7z(R, t) = p(R, R, t) - y -^ /(p, R, t). (67) 

Similarly, the current operator: {l/2mi) (t/'^V?/' — (VV'^)'0) has for expecta- 
tion value: 

j(R,t) = -^(V,-V,)p(y,x,i)||y_,|^o= /7^-/(p,R,i).(68) 
2mi J {2-Ky m 

These results are indeed those one would obtain in a classical description with 
/(p, R, t) the probability density to find a particle with momentum p at point 
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R and time t. Note however that while / is real, due to the hermiticity of p, 
it is not always positive as a truly classical distribution function would be. Of 
course / contains the same quantum information as p, and it does not make 
sense to specify quantum mechanically both the position and the momentum. 
However, / behaves as a classical distribution function in the calculation of 
one-body observables for which the typical momenta p that are involved in 
the integration are large in comparison with the scale 1/A characterizing the 
range of spatial variations of /, i.e. pX^ 1. 

By using the equations of motion for the field operators, i!Z'(r, t) = [H, 'F], 
where H is the single particle Hamiltonian, one obtains easily the following 
equation of motion for the density matrix 

idtp=[H,p]. (69) 

In fact we shall need the Wigner transform of this equation in cases where 
the gradients with respect to R are small compared to the typical values of 
p. Under such conditions, the equation of motion reduces to 

^^f + VpH-VRf-VRH-Vpf = 0. (70) 

where we have kept only the leading terms in an expansion in Vi^. For par- 
ticles interacting with gauge potentials A^{X), the Wigner transform of the 
single particle Hamiltonian in (|7y) takes the form: 

2 2 

H(R, p, f) = |- - -A • p + -A2(R, t) + eAo(R, t). (71) 

Zm m m 

Assuming that the field is weak and neglecting the term in A^ , one can write 
@l in the form: 

5t/ + V • VrJ + e(E + V a B) • Vp/ + ^{p^d.A')^^! = 0, (72) 

where we have set v = (p — eA)/m. This equation is almost the Vlasov 
equation (^5|): it differs from it by the last term which is not gauge invariant. 
The presence of such a term, and the related gauge dependence of the Wigner 
function, obscure the physical interpretation. It is then convenient to define 
a gauge invariant density matrix: 

/5(r,r',i) = (7^t(i.',i)^(r,i))C/(r,r',i), (73) 

where (s = r — r') 

t/,.O^exp(-./;dz^A,z..„)»expH.s^A(R„ ,74) 

and the integral is along an arbitrary path going from r' to r. Actually, in 
the last step we have used an approximation which amounts to chose for this 
path the straight line between r' to r; furthermore, we have assumed that the 



Theory of the quark-gluon plasma 21 

gauge potential does not vary significantly between r' to r. (Typically, p{r, r') 
is peaked at s = and drops to zero when s ^ At where At is the thermal 
wavelength of the particles. What we assume is that over a distance of order 
At the gauge potential remains approximately constant.) Note that in the 
calculation of the current, only the limit s — *■ is required, and that is given 



correctly by (74) (see also ( |75[ ) below). With the approximate expression ( [74[ ) 
the Wigner transform of (O) is simply /(R, k) = /(R, k + eA). By making 
the substitution /(R, p) = /(R, p — eA) in ([72|), one verifies that the non 
covariant term cancels out and that the covariant Wigner function / obeys 
indeed Vlasov's equation. 

In the presence of a gauge field, the previous definition ( |68| ) of the current 
suffers from the lack of gauge covariance. It is however easy to construct a 
gauge invariant expression for the current operator. 



'4>\-V - eA)i) - ({-V + eA)iljA ijA , (75) 



2m 



whose expectation value in terms of the Wigner transforms reads: 

The last expression involving the covariant Wigner function makes it clear 
that j(R, i) is gauge invariant, as it should. The momentum variable of the 
gauge covariant Wigner transform is often referred to as the kinetic momen- 
tum. It is directly related to the velocity of the particles: k = mv = p — eA. 
As for p, the argument of the non-covariant Wigner function, it is related to 
the gradient operator and is often referred to as the canonical momentum. 

In order to understand the structure of the equations that we shall obtain 
for the QCD plasma, it is finally instructive to consider the case where the 
particles possess internal degrees of freedom (such spin, isospin, or colour). 
The density matrix is then a matrix in internal space. As a specific example, 
consider a system of spin 1/2 fermions. The Wigner distribution reads |37[ : 

/(p,R) = /o(p,R) + /a(p,R)a„ (77) 

where the <Ja are the Pauli matrices, and the /q are three independent dis- 
tributions which describe the excitations of the system in the various spin 
channels; together they form a vector that we can interpret as a local spin 
density, f = (l/2)Tr(/cr). When the system is in a magnetic field with Hamil- 
tonian H = — /zq <t • B the equation of motion for f acquires a new component, 
9tf = 2/zoB A f , which accounts for the spin precession in the magnetic field. 
In the linear approximation this precession may be viewed as an extra time 
dependence of the distribution function along the characteristics: 

^ = ^ + V • Vtj + 2/ioB A . (78) 
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It is important to realize that all the differential operators above and in the 
Vlasov equation apply to the arguments of distribution functions, and not to 
the coordinates of the actual particles. Note however that equations similar to 
the ones presented here can be obtained for classical spinning particles. When 
the angular momentum of such particles is large, it can indeed be treated as a 
classical degree of freedom, and the corresponding equations of motion have 
been written by Wong |3^. After replacing spin by colour, these equations 
have been used by Heinz |[39| , [40|] in order to write down transport equations 
for classical coloured particles. By implementing the relevant kinematic ap- 
proximations one then recovers [Q the non-Abelian Vlasov equations to be 
derived below, i.e., (|7|) and i^. (See also |||,|^ for related work.) 



5.3 QCD Kinetic equations and hard thermal loops 

We are now ready to present the equations that are obtained for the QCD 
plasma. These are equations for generalized one-body density matrices de- 
scribing the long wavelength collective motions of colour particles (quarks 
and gluons), and possible excitations involving oscillations of fermionic de- 
grees of freedom. They look formally as the Vlasov equation, the main ones 
being (§11 : 

[vD,,Sn±{k,x)]^T9^-'E{x)^, (79) 

[vD,,SN{k,x)]=~g^-E{x)'^, (80) 

(v ■ I?,)/(k, x) = -igCf {Nk + Uk) i ^(x). (81) 

In these equations, v^ — (l,v), v = k/fc, '^{x) is an average (relativistic) 
fermionic field, and (5n±, SN and yl are gauge-covariant Wigner functions for 
the hard particles. The first two Wigner functions are those of the density 
matrices of the quarks and the gluons, respectively; the last one is that of 
a more exotic density matrix which mixes bosons and fermions degrees of 
freedom, A ^ {ipA). The right hand sides of the equations specify the quan- 
tum numbers of the excitations that they are describing: gluon for the first 
two, quark for the last one. One of the major difference between the QCD 
equations above and the linear Vlasov equation for QED is the presence of 
covariant derivatives in the left hand sides of the equations. These play a 
role similar to that of the magnetic field in ( [78|) for the distribution functions 
of particles with spin. (Note that the equation for / holds for QED, with a 
covariant derivative there as well.) 

The equations (|79[)-(|8l|) have a number of interesting properties which 
are reviewed in |l||. In particular, they are covariant under local gauge trans- 
formations of the classical fields, and independent of the gauge-fixing in the 
underlying quantum theory. 
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By solving these equations, one can express the induced sources as func- 
tional of the background fields. To be specific, consider the colour current: 

j^{x) ^ 2.9 1 ^ v^ Tr (T'^SNik, x)) , (82) 

where SN is the gluon density matrix. Quite generally, the induced colour 
current may be expanded in powers of A^ , thus generating the one-particle 
irreducible amplitudes of the gauge fields |^ : 

j;i = i7^X + i^;.V^^c+- (83) 

Here, i7°[j = 6'^'^n^i, is the polarization tensor, and the other terms repre- 
sent vertex corrections. These amplitudes are "hard thermal loops" (HTL) 
[ p0pl| , |3^ , p5p^ which define the effective theory for the soft fields at the scale 
gT . It is worth noticing that the kinetic equations isolate directly these hard 
thermal loops, in a gauge invariant manner, without further approximations. 
The gluon density matrix can be parametrized as in ( pq ) : 

SNabiK^) = -gWab{x,^){ANk/dk), (84) 

where Nk = l/(e'^'' — 1) is the Bose-Einstein thermal distribution, and 
W{x, v) = Wa{x, v)r" is a colour matrix in the adjoint representation which 
depends upon the velocity v — k/fc (a unit vector), but not upon the mag- 
nitude k — |k| of the momentum. Then the colour current takes the form: 

jLM-^IJ^^^W{x,^) (85) 

with m?jj ~ g^T^ ■ A similar representation holds for the quark density matri- 
ces 5n± (k, x) . The kinetic equations for 5Nab and 5n± can then be written 
as an equation for Wa{x,w): 

(vi^,)"^VFb(x,v) =v-E'^(a;). (86) 

They differ from the corresponding Abelian equation ( |60| ) merely by the 
replacement of the ordinary derivative dx ~ gT by the covariant one D^ = 
dx + igA. Accordingly, the soft gluon polarization tensor derived from (p5|)- 
(p6|), i.e., the "hard thermal loop"_/7^^, is formally identical to the photon 
polarization tensor obtained from (^0^ and given by (|6^) |^,^ . The reason 
for the existence of an infinite number of hard thermal loops in QCD is the 
presence of the covariant derivative in the left hand side of (p6|). A similar 
observation can be made by writing the induced electromagnetic current in 
the form: 

jU^)=mll^v'^Jd^y{x\-^\y)v-E{y) 

= Jd*ya^'^x,y)E^iy). (87) 



24 Jean-Paul Blaizot 

This expression, which is easily obtained from the expression (m) of Sf, 
defines the conductivity tensor a^'^ . The generahzation of this expression to 
QCD amounts essentiahy to replacing the ordinary derivative by a covariant 
one. 



6 Collective phenomena in the quark-gluon plasma 

At the classical level, the effective theory at the scale gT is summarized by 
the generalized Yang-Mills equations 



in V ■ D 



D^F'^^ ^nil I ^^ E^ ^ nf^A-^ + - rXAlA<^, + ... (88) 



where the induced current in the right hand side describes the polarization 
of the hard particles by the soft colour fields A(^. In this equation, rhn ~ gT 
is the Debye mass, E\ is the soft electric field, v^^ = (1, v), and the angular 
integral J AQ runs over the orientations of the unit vector v. The induced 
current is non-local and gauge symmetry, which forces the presence of the 



covariant derivative D^^ = d'^ + igA'^ in the denominator of (88), makes it 
also non-linear. 

Similarly, the soft fermionic fields obey the following generalized Dirac 
equation l^^l (with M ^ gT and ;*^ = 7pu'') : 



. , dn ^ 



,p^ -M^J — -^^ ^^Sr^ + T'^A^^i, + ... (89) 

These equations allow the description of a variety of collective phenomena. 
We discuss briefly here some of them ( collective modes, Debye screening and 
Landau damping). More details can be found in the lecture by A. Rebhan. 
See also ilT 



6.1 Collective modes 

The collective plasma waves are propagating solutions to (pq) or (p9|) . We re- 
strict ourselves in this subsection to the weak field limit where these equations 
become linear and essentially Abelian. 

The solutions can then be analyzed with the help of the propagator. We 
consider here the gluon propagator *G^i,, in Coulomb's gauge where it has the 
following non-trivial components, corresponding to longitudinal (or electric) 
and transverse (or magnetic) degrees of freedom: 

*Gooito,p) = *Al{uj,p), *G,,(a;,p) = (5,, -p,;p,)*Z\t(c^,p), (90) 

where: 

*Al{uj,p) = ~\ r, *ATiLJ,p) = — T-^n V' (91) 

p^ + nL[u),p) Lu-^ -p^ - nT[uj,p) 
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and the electric {Ul) and magnetic (^7^) polarization functions are defined 
as: 

nL{co,p) = -noo{co,p) , nT{oj,p) = \ {6^^ - p'fP)n,j[Lo, p) . (92) 

Explicit expressions for these functions can be found in ||l|. 





(a) (b) 

Fig. 6. Dispersion relation for soft excitations in the linear regime: (a) soft fermions; 
(b) soft gluons (or linear plasma waves), with the upper (lower) branch correspond- 
ing to transverse (longitudinal) polarization. 



The dispersion relations for the modes are obtained from the poles of the 
propagators, that is. 



P 



iTi(c^L,p) = 0, 



2 

LOrp 



P 



nT{ujT,p), 



(93) 



for longitudinal and transverse excitations, respectively. The solutions to 
these equations, ujl{p) and Ci;r(p), are displayed in Fig. o.b. The longitudinal 
mode is the analog of the familiar plasma oscillation. It corresponds to a col- 
lective oscillation of the hard particles, and disappears when p ^ gT . Both 
dispersion relations are time-like (lol^t{p) > p), and show a gap at zero mo- 
mentum (the same for transverse and longitudinal modes since, when p — j 0, 
we recover isotropy). With increasing momentum, the transverse branch be- 
comes that of a relativistic particle with an effective mass ttjoo = mo/V^ 
(commonly referred to as the "asymptotic mass"). Although, strictly speak- 
ing, the HTL approximation does not apply at hard momenta, the above 
dispersion relation ujt{p) remains nevertheless correct for p ^ T where it 
coincides with the light-cone limit of the full one-loop result |H : 



i7^-'°°P(co2=p2) ^ !!^. 



(94) 
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The dispersion relations of soft fermionic excitations exhibit also collective 
feature with a characteristic splitting at low momenta (see Fig. 0.a). We shall 
not discuss here this interesting phenomenon (see B and references therein). 

We note finally that particular solutions of the non-linear equations ( p8| ) 
have also been found, in |43,Eq,H. These solutions describe non-linear plane- 
waves propagating through the plasma, and represent truly non-Abelian col- 
lective excitations. 

6.2 Debye screening 

The screening of a static chromoelectric field by the plasma constituents is 
the natural non-Abelian generalization of the Debye screening, a familiar 
phenomenon in classical plasma physics |33[ . In coordinate space, screening 
reduces the range of the gauge interactions. In momentum space, it con- 
tributes to regulate the infrared behaviour of the various n-point functions. 
Screening properties can be inferred from an analysis of the effective pho- 
ton (or gluon) propagators (|9l|) in the static limit uj —> 0. We have: 



7^D> 


7Tt(0,p) = 0, 




-1 


*^t(0,p) = 


1 


p2 + ml, 


- ^ 



nLiO,p)^ml, 7Tt(0,p) = 0, (95) 

and therefore: 

*Al{0,p) = 2^^ 2 , *Z\t(0,p)--, (96) 

p-^ + mjj p'^ 

which clearly shows that the Debye mass mB acts as an infrared cut-off ~ gT 
in the electric sector, while there is no such cut-off in the magnetic sector. 

6.3 Landau damping 

For time-dependent fields, there exists a different screening mechanism asso- 
ciated to the energy transfer to the plasma constituents. In Abelian plasmas, 
this mechanism is known as Landau damping |33[ . The mechanical work done 
by a longwavelength electromagnetic field acting on the charged particles 
leads to an energy transfer |33| : 



dEM/(t) 
dt 



= yd3xE(i, x).j(i, x), (97) 



where f{p) — n^^{p)A^{p) is the induced current. One can then show that 
the average energy loss is related to the imaginary part of the retarded po- 
larization tensor. From (p3h we get: 

lmn^''{uj,p) ^ -nmlu; f ^ v^'v''6{iJ-v■p). (98) 

The (5-function in (BSf) shows that the particles which absorb energy are those 
moving in phase with the field (i.e., the particles whose velocity component 
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along p is equal to the field phase velocity: v • p = io/p). Since in ultrarela- 
tivistic plasmas v is a unit vector, only space-like (|w| < p) fields are damped 
in this way. 

To see how this mechanism leads to screening, consider the effective pho- 
ton (or gluon) propagator (p^), and focus on the magnetic propagator. For 
small but non-vanishing frequencies the corresponding polarization function 
iTr(a;,p) is dominated by its imaginary part: 

nriuj <^p) = -t^ml-+ Oico^p') , (99) 

4 p 

and therefore 

*ATiu; « p) ~ -, ., \^ , , . (100) 

Thus ImiTT(p) acts as a frequency-dependent IR cutoff at momenta p ^ 
{iVTrij^y/'^. That is, as long as the frequency uj is different from zero, the soft 
momenta are dynamically screened by Landau damping p7[|. 



7 The entropy of the quark-gluon plasma 

We come now to the last part of these lectures which will be mainly devoted 
to an introduction to the recent progress made in the calculation of the 
entropy of the quark-gluon plasma. We first comment on various aspects of 
perturbation theory and show that it is not appropriate for calculating the 
thermodynamics of the quark-gluon plasma, even a high temperature where 
the coupling is weak. The main source of difficulties is that the contributions 
of the collective modes, for which we have constructed an effective theory 
in the previous sections, are non perturbative and cannot be expanded in 
powers of the coupling constant. We then show that these contributions can 
be included by using self-consistent approximations familiar in many-body 
physics. These are best formulated for the entropy of the plasma, for which 
we obtain a simple approximation which provides an accurate description of 
lattice gauge calculations. 

7.1 Results from perturbation theory 

The free energy has been calculated up to order g^ , including the contribution 
of fermions [Bq|. However, since our purpose here is mostly pedagogical, we 
shall limit our discussion to the gluon contribution at order g"*, in an SU(N) 
gauge theory. The pressure P = —F/V can then be written: 

P = Po[l + a2.g' + as.?' + (aiWT) + a'^ ln.g) g^ + Oig^)] , (101) 
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(102) 

where C is Riemann's zeta function, and /i the renormaUzation scale. 

The first term in the expansion is Pq, the pressure of an ideal gas of 
gluons: 

Po = (N^ - 1) T^j!. (103) 

The next term, of order (7^, is a genuine perturbative correction, and so is the 
term of order g*. The contributions of order g^ can be interpreted as a contri- 
bution of the collective modes to the pressure, and the odd power reflects the 
fact that the calculation of this contribution requires resummations. Similar 
resummations are responsible for the term in g'^lng. 





2 S / 2 da ^ 
M 2 + f^ A 2 


d 
da 


In th 
satis 


is equation, a(/i) e 
ies the equation: 


^9' 



We note that some of the coefficients in (102) depend on the renormaliza- 
tion scale /z. However, the pressure itself should not depend on /i. It obeys a 
renormalization group equation: 

P = 0. (104) 

(7^(/x)/47r is the running coupling constant which 
e equation: 

A*'^ - /3(a) = -/32a' -/33a', (105) 

with 

117V IIN"^ 

P^ = Tr^ ^3 = ^. (106) 

127r 247r^ 

One can then show that, indeed, P is independent of /i: the explicit /x de- 
pendence of the coefficients cancels with that of the running coupling. Look 
indeed at the following combination of terms coming from the contributions 
of 025^ and the ^, dependent part of 04^^: 

N \ ^ , N 2 22,_^ ^l 



4^1"+4^"yi^4^^ (^«^) 

By taking the derivative of this expression with respect to ^^ one gets: 

n d „ n da A^ 9 11 , . , , /,„„N 

^ -— 5- {} = ^ -— r H a^ h higher order terms. (108) 

d/z^ d^^ 47r 3 
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By using the leading order expression of the /3- function given in (JlOq), one 
then obtains, as announced: 

Note however, that the pressure is only formally independent of /i at order 
g^, in the sense that its derivative with respect to /i involves terms of order g^ 



at least. But the approximate expression (101) for P does depend on fi. As in 
all perturbative calculations, one is then led to look for the best value of /z, i.e. 
the one which minimizes the higher order corrections. In the present context, 
a "natural choice" is to fix /i = 2ttT, where 2ttT is the scale provided by the 
basic Matsubara frequency. This choice makes the running coupling decrease 
with increasing temperature, and leads in particular to the expectation that 
the quark-gluon plasma becomes perturbative at very l arge temperature. 



By cal culat ing explicitly the various coefficients in ( |102[ ) for A^ = 3, one 



can write (|101[ ) thus: 

P^Po[l- 0.095.9^ + 0.12.9^- 



(0.09 In^ 



0.007- O.OlSln (2^)) g^ + 0(g^)| . (110) 



Then, if for example one fixes fi — 2t:T and chooses a large temperature such 
that a(27rr) == 0.1, one gets g = 1.12, and 

P = Po [1-0.12 + 0.17 + 0.004], (HI) 

which shows no sign of convergence, with the term of order g^ larger than the 
term of order g^. Furthermore, if one analyzes the dependence of P on the 
renormalization scale, on finds large variations as /i runs within the interval 
■kT < n< AttT. 

Attempts have been made to extract information from the first terms of 



this series using Fade approximants [£3 54 or Borel summation techniques 
[ p5|j5q| . The resulting expression of the pressure becomes indeed a smooth 
function of the coupling, better behaved than the polynomial approximation 



(101). These techniques however, which are in some situations very powerful, 
provide little physical insight, and we shall not discuss them further here. 

The behavior of perturbation theory does not improve as one takes into 
account the higher order terms that one can calculate (namely orders g^ 
and g^). Furthermore, at order g^, as we have already mentioned, perturba- 
tion theory becomes inapplicable because of infrared divergences. It has been 
shown in [ p^ , poH5l] ] how, in principle, an effective theory could be constructed 
to overcome this particular problem by marrying analytical techniques (to de- 
termine the coefficients of the effective theory) and numerical ones (to solve 
the non perturbative 3-dimensional effective theory). The resulting effective 
theory is a 3-dimensional theory of static fields, with Lagrangian: 

^eff = \{Ft,r + liD.A^of + lrnl{A^,r + A(AS)^ + 6C, (112) 
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with Di = di—igvTAi. This strategy has been apphed recently to the calcula- 
tion of the free energy of the quark-gluon plasma a high temperature |53] . The 
slow convergence of the pressure towards the ideal gas value, that is seen in 
lattice calculations above Tc, is well reproduced. It is worth-emphasizing that 
this technique of dimensional reduction puts a special weight on the static 
sector (it singles out the contributions of the zero Matsubara frequency). 
However, as we shall see, it may be advantageous to keep, even in the calcu- 
lation of equilibrium thermodynamic properties, the full spectral information 
that one has about the plasma excitations. 

There are indeed indications that lattice data are well accounted for by 
simple phenomenological models of weakly interacting quasiparticles |57| , p8| . 
In the case of the scalar field, the dominant effect of the interactions is to 
give a mass to the excitations. An indeed a perturbative expansion in terms 
of screened propagators (that is keeping the screening mass '^ gT as a pa- 
rameter, not considered as a term of order g entering the expansion) has been 
shown to be quite stable with good convergence properties ^. In the case 
of gauge theory, the effect of the interactions is more complicated than just 
generating a mass. But we know how to determine the dominant corrections 
to the self- energies. When the momenta are soft, these are given by the hard 
thermal loops discussed above. By adding these corrections to the tree level 
Lagrangian, and subtracting them from the interaction part, one generated 
the so-called hard thermal loop perturbation theory Q . The resulting per- 
turbative expansion is made complicated however by the non local nature of 
the hard thermal loop action, and by the necessity of introducing tempera- 
ture dependent counter terms. At the expense of some extra formalism, some 
of these difficulties can be avoided. This is discussed now. 



7.2 Skeleton expansion for thermodynamic potential and entropy 

In this section we recall the formalism of propagator renormalization that al- 
low systematic rearrangements of the perturbative expansion while avoiding 
double-counting. We shall see in particular how self-consistent approxima- 
tions can be used to obtain a simple expression for the entropy which isolates 
the contribution of the elementary excitations as a leading contribution. For 
pedagogical purposes, we shall mainly consider in these lectures the example 
of the scalar field. 

The thermodynamic potential i7 = —PV of the scalar field can be written 
as the following functional of the full propagator D [pllp2l : 



f3f2[D] = -\ogZ = -Tr \ogD-^ - -TrnD + <P[D], (113) 

where Tr denotes the trace in configuration space, /3 = l/T, 77 is the self- 
energy related to D by Dyson's equation (Dq denotes the bare propagator): 

D-^^Dq^+H, (114) 
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and ^[D] is the sum of the 2-particle-hTeducible "skeleton" diagrams 



<P[D]= 1/12^— ^ +1/8 r'V^ + 1/48 r^+--- (115) 



The essential property of the functional f2[D] is to be stationary under 
variations of D (at fixed Dq) around the physical propagator. The physical 
pressure is then obtained as the value of f2[D] at its extremum. The station- 
arity condition, 

Sn[D]/5D = 0, (116) 

implies the following relation 

S^[D]/SD ^ iyj, (117) 



which, together with (114), defines the physical propagator and self-energy 



in a self-consistent way. The equation (117) expresses the fact that the skele- 



ton diagrams contributing to U are obtained by opening up one line of a 
two-particle-irreducible skeleton. Note that while the diagrams of the bare 
perturbation theory, i.e., those involving bare propagators, are counted once 
and only once in the expression of 11 given above, the diagrams of bare per- 
turbation theory contributing to the t herm odynamic potential are counted 
several times in <P. The extra terms in (|ll3| ) precisely correct for this double- 
counting. 

Self-consistent (or variational) a ppro ximations, i.e., approximations which 



preserve the stationarity property (116), are obtained by selecting a class of 



skeletons in 'P[D] and calculating U from (117). Such approximations are 
commonly called "^-derivable" ||62|. 



The traces over configuration space in (113) involve integration over imag- 
inary time and over spatial coordinates. Alternatively, these can be turned 
into summations over Matsubara frequencies and integrations over spatial 
momenta: 

At I d?x~^ fiV ([dkl (118) 

where V is the spatial volume, k^ = (ici;„,k) and w„ — mrT, with n even 
(odd) for bosonic (fermionic) fields (the fermions will be discussed later). 
We have introduced a condensed notation for the the measure of the loop 
integrals (i.e., the sum over the Matsubara frequencies ijJ„ and the integral 
over the spatial momentum k): 



n.even 



Strictly speaking, the sum-integrals in equations like (113) contain ultraviolet 



divergences, which requires regularization (e.g., by dimensional continuation). 
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Since, however, most of the forthcoming calculations will be free of ultraviolet 
problems, we do not need to specify here the UV regulator (see however 



Sect. 7.3 for explicit calculations). 

For the purpose of developing approximations for the entropy it is con- 
venient to perform the summations over the Matsubara frequencies. One ob- 
tains then integrals over real frequencies involving discontinuities of propaga- 
tors or self-energies which have a direct physical significance. Using standard 
contour integration techniques, one gets: 

njV = Ittt^ "(^) (lmlog(--w2 + k^ + n) - IninD) + T<P\D\IV 

(120) 

where n(w) = l/(e'''^ - 1). 

The analytic propagator D{uj, k) can be expressed in terms of the spectral 
function: 

D[u.^k)= r ^Pi^iR, (121) 

J-oo 27r ko-uj 

and we define, for uj real, 

ImD(cj,k)=ImD(cj + ie,k) = ^^^. (122) 

The imaginary parts of other quantities are defined similarly. 
We are now in the position to calculate the entropy density: 

S = -d{f2/V)/dT. (123) 



The thermodynamic potential, as given by ( 120 ) depends on the tempera- 



ture through the statistical factors n{uj) and the spec tral function p, which 



is determined entirely by the self-energy. Because of ( 116 ) the temperature 
derivative of the spectral density in the dressed propagator cancels out in the 
entropy density and one obtains [|63H64| : 



5--y(^^ImlogD-(.,k) 
df^k dn{uj) 



with 
S' 



(27r)4 dT 
d{T$) f d*k dn{uj 



Imn(w, k)ReD(w, k) + S' (124) 



dT 



RenimD. (125) 



D J (27r)4 dT 
For the two-loop skeletons, we have: 

S' = 0. (126) 
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Loosely speaking, the first two terms in (124) represent essentially the entropy 
of "independent quasiparticles" , while S' accounts for a residual interaction 
among these quasiparticles [ |64| . 

The vanishing of iS' holds whether the propagator are the self-consistent 



propagators or not. That is, only the relation (117) is used in the proof 
which does not require D to satisfy the self-consistent Dyson equation (114). 
A general analysis of the contributions to S' and their physical interpretation 
can be found in |p5|. 



We emphasize now a few attractive features of the formula (124) with 
S' = 0, which makes the entropy a privileged quantity to study the thermo- 
dynamics of ultrarelativistic plasmas. We note first that the formula for S 
at 2-loop order involves the self-energy only at 1-loop order. Besides this im- 
portant simplification, this formula for S, in contrast to the pressure, has the 
advantage of manifest ultra-violet finiteness, since dn/dT vanishes exponen- 
tially for both Lo — > ±cxD. Also, any multiplicative renormalization D -^ ZD, 



n ^ Z^^n with real Z drops out from (124). Finally, the entropy has a more 
direct quasiparticle interpretation than the pressure. This will be illustrated 
explicitly in the simple model of the next subsection. 

7.3 A simple model 

In this section we shall present the self-consistent solution for the (A/4!)(/)^ 
theory, keeping in (p only the two-loop skeleton. Anticipating the fact that 
the fully dressed propagator will be that of a massive particle, we write the 
spectral function as p(fco, k) = 27r e(A;o) 5{kQ — k^ — m^), and consider m as a 
variational parameter. The thermodynamic potential ( |113| ), or equivalently 
the pressure, becomes then a simple function of to. By Dyson's equation, the 
self-energy is simply 77 — m? . We set: 

/(m) ^ \ f[dk] D{k) = i /[dfc] ,^.1^ , . (127) 

2 ./ 2 ./ uji, + k"^ + TO-^ 



Then the pressure P = —il/V can be written as: 



1 f d^fc 1 f d^fc , , o^ , 9^, . Ao o, 

2y(^^^-^W(^^°^^'-^^^""^'^™^+f 

(128) 

where ef = k'^ + m?. By demanding that P be stationary with respect to 
m one obtains the self-consistency condition which takes here the form of a 
"gap equation" : 

to2 ^ Ao/(m). (129) 



The pressure in the two-loop <?-derivable approximation, as given by (127) 



(129), is formally the same as the pressure per scalar degree of freedom 



in the (massless) A'^-component model with the interaction term written as 



34 Jean-Paul Blaizot 

j^r^i'^/ 4:l){4'i(t>i)^ in the limit A^ ^ oo [p6|. From the experience with this lat- 
ter model, we know that ( |l27| )-( p^ ) admit an exact, rcnormalizable solution 
which we recall now. 

At this stage, we need to specify some properties of the loop integral 
I{m) which we can write as the sum of a vacuum piece Io{m) and a finite 
temperature piece Irim) such that, at fixed m, Irim) ^- as T — > 0. We 
use dimensional regularization to control the ultraviolet divergences present 
in Iq, which implies /o(0) = 0. Explicitly one has: 

^^/(m) = -^ (- + log 4 + l) + Irim) + 0(e), (130) 



327r2 V e m'^ 



with 






and Ek = {k^ + m.^)^/^. In (13C), /i is the scale of dimensional regulariza- 
tion, introduced, as usual, by rewriting the bare coupling Aq as /i'^Ao, with 
dimensionless Aq; furthermore, e = 4 — n, with n the number of space-time 
dimensions, and p? = 47re~'''/i^. 

We use the modified minimal subtraction scheme (MS) and define a di- 
mensionless renormalized coupling A by: 

a = a;;;^ + i6^- ^^^^^ 

When expressed in terms of the renormalized coupling, the gap equation 
becomes free of ultraviolet divergences. It reads: 



A f d^fc n{ek) Am^ 



m 



2 



"^' - o t:^ ^r^ + ^^^^ logTT - 1 , (133) 



2 J (27r)3 Ek 327r2 V M^ 

The renormalized coupling constant satisfies 
dA A2 



dlog^ IGtt^ 



(134) 



whic h ens ures that the solution m of (133) is independent of /i. The expres 



sion (134) coincides with the exact /3- function in the large- A^ limit, but gives 
only one third of the lowest-order perturbative /3-function for A^ = 1. This 
is no actual fault since the running of the coupling affects the thermody- 
namic potential only at order A^ which is beyond the perturbative accuracy 
of the 2-loop ^-derivable approximation. In order to see the correct one-loop 
/3-function at finite N, the approximation for <P would have to be pushed to 
3-loop order. 



Note also that, in the present approximation, the renormalization (132) of 



the coupling constant is sufficient to make the pressure (128) finite. Indeed 
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in dimensional regularization the sum of the zero point energies 6^/2 in (128) 
reads: 



(135) 



"V 


(27r)"-i 2 


m4 


- + log ^ -f 


647r2 


so that 








"■y 


(27r)"-i 2 


7T2 _ 

2Ao 


4 4 


2A 647r2 



+ K +0(e), 



l°g4 + ?)+0W (136) 



is indeed UV finite as n -^ 4. After also using the gap equation (|133| ), one 
obtains the /l-independent result 



P 



-T 



d^fc 

^2^ 



log(l 



-fisk 



) + -7rlT{m) + 



1287r^ 



(137) 



We now compute the entropy according to (124). Since Imll — and 
Rell = ni^, we have simply: 



S 



d^k dn{Ld) 



Imlog(k^ -u"^ + T0? 



(27r)4 dT 
Using 

Imlog(k2 - w^ + m^) ^ -Tre{Lj)9{Lj^ ~ el), 
and the identity, 

dn{uj) da{uj) 



dT 



duj 



cr{ijj) = -71 log 71+ (1 + rt)log(l + n), 



one can rewrite (138) in the form (with n^ = n{ek))'- 
5 = 



<Pk 

(27r) 



— ^|(l + nfe)log(l + 7ifc) - nklogukj. 



(138) 



(139) 



(140) 



(141) 



This formula shows that, in the present approximation, the entropy of the 
interacting scalar gas is formally identical to the entropy of an ideal gas of 
massive bosons, with mass m. 

It is instructive to observe that such a simple interpretation does not hold 
for the pressure. The pressure of an ideal gas of massive bosons is given by: 



P^°\m) 



(Pk 

(2^.e. 



dw ( n{ijj) + 



[ (Pk ( 



1 



"'-)+l}, 



(142) 



which differs indeed from (128) by the term rn^/X which corrects for the 
double-counting of the interactions included in the thermal mass. 
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7.4 Comparison with thermal perturbation theory 

In view of the subsequent application to QCD, wliere a fully self-consistent de- 
termination of the gluonic self-energy seems prohibitively difficult , we shall be 
led to consider approximations to the gap equation. These will be constructed 
such that they reproduce (but eventually transcend) the perturbative results 
up to and including order A'^" or g^, which is the maximum perturbative 
accuracy allowed by the approximation S' = 0. 

In view of this it is important to understand the perturbative content of 
the self-consistent approximations for m^, P and S. In this section we shall 
demonstrate that, when expanded in powers of the coupling constant, these 
approximations reproduce the correct perturbative results up to order A'^" 
[ [ll| . This will also elucidate how perturbation theory gets reorganized by the 
use of the skeleton representation together with the stationarity principle. 

For the scalar theory with only (A/4!) (j)^ self-interactions, we write|J A = 
24(7^, and compute the corresponding self-energy 77 = rn^ by solving the 
gap equation ( |l33| ) in an expansion in powers of 5, up to order g^. Since we 
anticipate m to be of order gT, we can ignore the second term ex Xm^ ~ g in 
the r.h.s. of ( |133[ ), and perform a high-temperature expansion of the integral 
Iri'm) in the first term (cf. ( |l3l| )) up to terms linear in m. This gives the 
following, approximate, gap equation: 

m^ ~ g'^T^ g^Tm. (143) 

TT 

The first term in the r.h.s. arises as 

24g2/T(0) = 12ff^/(0 ^ = 9'T' ^ rn\ (144) 

This is also the leading-order result for rn^, commonly dubbed the "hard 
thermal loop" . 



The second term, linear in m, in (143) comes from 
19 ^ / "^^^ fn{sk) n{ky 



(2^)3 V Sk 



(145) 



d^k ( 1 1 ^^ _ 3g2 

(2^ \k'^+m? " fc2 



iVr Ij-r, ( i:^^3 - 1:2 ) = -^^T, (146) 



where we have used the fact that the momentum integral is saturated by soft 
momenta k ~ gT, so that to the order of interest n{ek) — T /eu (and similarly 

^ This normalization for g is cliosen in vie w of t lie subsequent extension to QCD 



since it makes tiie scalar thermal mass in ( 144 ) equal to the leading-order Debye 
mass in pure-glue QCD. 
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n{k) ~ T/k). This provides the next-to-leading order (NLO) correction to the 
thermal mass 



5m 



ig' 



mT 



3t^2 



■g'T 



(147) 



Thus, to order g'^, one has m? = m^ + Sm^. In standard perturbation 
theory [0,0, the first term arises as the one- loop tadpole diagram evaluated 
with a bare massless propagator, while the second term comes from the same 
diagram where the internal line is soft and dressed by the HTL, that is 

Consider similarly the p crtur bativ e esti mates for the pressure and entropy, 
as obtained by evaluating ( 128 ) and (|l4l| ) with the pcrturbative self-energy 



n — 771 ~ m + Sm , an d fur ther expanding in powers of g, to order g . The 
renormalized version of (128) yields, to this order (recall that m '^ gT and 



2rp4 



P 



n^T 



m^T^ nv'T 



m 
2A 



90 24 ■ 127r ' ' 2A ' ^^^^' 

The first terms before the dots represent the pressure of massive bosons, 
i.e. (142) expanded up to third order in powers of m/T. From (148), it can 
be easily verified that the above pcrturbative solution for m? ensures the 
stationarity of P up to order g"^, as it should. Indeed, if we denote 



P2(m) = - 



24 



m 
2A 



then the following identities hold: 



dm 



= 0, 



dm 



Psim) 



dm 



m-^T 



(149) 



= 0. 



This shows that the NLO mass correction 6m,^ 
as 



(150) 
g^T^ can be also obtained 



3771 



(dPs/dm) 



{d^P2/dm^) 



3ff .2 

m . 

n 



(151) 



in agreement with (|l47| ). Moreover, P2 = P2{m) = —g'^T'^/48 and P3 



P^iiTi) = 7?7,^T/127r are indeed the correct pcrturbative corrections to the 
pressure, to orders g^ and g^, respectively ||ll|. In fact, the pressure to this 
order can be written as: 



P 



90 
90 



24 ^ TT ^' 



m-^T 
Utt 



Ir"-^)' 



777 
48 



T^ + 



Utt 



0{g') 
(152) 



Note that the term of order g^ is only /laZ/of that one would obtain from (142) 
by replacing 777 by m. This is due to the mismatch between (142) and the 
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correct expression (128) for the pressure. In fact the net order g^ contribution 
to the pressure comes from <1> evaluated with bare propagators: the order g^ 
contributions in the other two terms mutuaUy cancel indeed. This is to be 
expected: there is a single diagram of order g^; this is a skeleton diagram, 
counted therefore once and only once in <?. Observe also that the terms of 
order g'^ originating from the terms rh? and m^ mutually cancel; that is, the 
NLO mass correction 5m drops out from the pressure up to order g'^ . This is 
no accident: the c ance llation results from the stationarity of P at order 5^, 



the first equation ( |150| ) 



Consider now the entropy density. The correct perturbative result up 
to order g '^ m ay be obtained directly by taking the total derivative of the 
pressure, ( |152| ) with respect to T. One then obtains: 

4 /^2y4 ^2^2 ^Sj^X 



We wish, however, to proceed differently, using (141), or equivalently, 



since dP/dm = when m is a solution of the gap equation, by writing: 

(154) 



^_dP 



This yields: 

5 = - + ]+0(m^T), (155) 

T V 90 48 487r y ^ ' '^' ^ ' 

which coincides as expected with the expression obtained by expanding the 
entropy (141) of massive bosons, up to order (m/T)^. If we now replace m 



by its leading order value m, the resulting approximation for S reproduces 
the perturbative effect of order ~ g^, but it underestimates the correction of 
order g^ by a factor of 4. This is corrected by changing m to rh + 5m with 



5m — —"igrnjI-K in the second order term of (155). Note that although it 
makes no difference to enforce the gap equation to order g^ in the pressure 
(because of the cancellation discussed above), there is no such cancellation 
in the entropy. 

7.5 Approximately self-consistent solutions 

As we have seen, the 2-loop ^-derivable approximation provides an expression 
for the entropy 5 as a functional of the self-energy 77 which has a simple quasi- 
particle interpretation and is manifestly ultraviolet finite for any (finite) 77. 



These attractive features of the formula (124) are independent of the specific 
form of the self-energy, and will be shown to hold in QCD as well. Of course, 
within this approximation, the self-energy is uniquely specified: by the sta- 
tionarity principle, this is given by the self-consistent solution to the one-loop 
gap equation. In the scalar 0^-model, it was easy to give the exact solution 
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to this equation. In QCD, however, it wiU turn out that a fully self-consistent 
solution is both prohibitively difficult (because of the non-locality of the gap 
equation), and not really desirable (because gauge symmetry implies rela- 
tions between the renormalization of the propagators and that of the vertices, 
and the present approximation deals only with propagator renormalization) . 
This leads us to consider approximately self- consistent resummations, which 
are obtained in two steps: (a) An approximation is con structed for the solu- 



tion n to the gap equation, and (b) the entropy ( 124 ) is evaluated exactly 
(i.e., numerically) with this approximate self-energy. While step (b) above is 
unambiguous and inherently non-perturbative, step (a), on the other hand, 
will be constrained primarily by the requirement of preserving the maximum 
possible perturbative accuracy, of order g^. In addition to that, we shall add 
the qualitative requirement that the approximation for 77, and the ensuing 
one for S, are well defined and physically meaningful for all the values of g 
of interest, and not only for small g — that is, for all the values of g where 
the fully self-consistent calculation makes sense a priori. Finally, in the case 
of QCD, relaxing the requirement of complete self-consistency allows us to 
construct gauge invariant approximations. 




(a) (b) (c) (d) 

Fig. 7. QCD skeletons contributing to <1> at 2-loop order. Wiggly, plain and dotted 
lines refer respectively to gluons, quarks and ghosts 



We shall now, in the rest of this lecture, outline the main steps that are 
involved in the implementation of these approximations in the case of QCD. 
Details can be found in the original publications [13 



14 



U 



At 2-loop order, the relevant skeletons are displayed in Fig. 7.5. By itself, 
the corresponding self-consistent truncation is not a gauge invariant approx- 
imation. Our strategy then will be to use gauge-invariant approximations 
to self-energies, in place of the self-consistent ones. These self-energies are 
then used to compute the entropy without further approximations. In com- 
plete analogy with the example of the scalar case that we have discussed in 
the previous section, these approximations are such that, when expanded in 
powers of the coupling the entropy is identical to that given by perturbation 
theory up to and including order g^. 

The approximate self-energies that we use are the hard thermal loops 
discussed above. Namely, for soft momenta ui, p ^ gT, we take TT^o/t ~ 
nuTL and Ssoft ~ ^htl, for gluons and quarks respectively. We shall also 
need an approximation valid for lu, p ^ T: Uhardi'^'^ ^ P^) and similarly for 
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^/vWw^ .AAMAAAA/OkAA/^ vAA/^AAA/'a^A/^ 



5nl 6nl snl, enl 



Y^ 



^^-* — » — u^ 

Fig. 8. Next to leading order contribution to SUt (top) and to SS (bottom) at hard 
momentum. Thick dashed and wiggly lines with a blob represent HTL-resummed 
longitudinal and transverse propagators, respectively 



S. It turns out that this is accurately given by the hard thermal loop, even 
though the momenta are not soft H]. All these approximations are gauge 



invariant. The corresponding diagrams are displayed in Figs. 71. 

We can then proceed exactly as in the scalar case. As a first approximation 
one may simply use the hard thermal loops 7T — IIhtl and E ^ Shtl 
at all momenta; we refer to the corresponding entropy as 5 = Shtl- The 
perturbative content of this approximation is schematically 0{g^) + , j 0{g^); 
that is, the approximation fully accounts or the order g^ , but reproduce only 
a quarter of the g^ order, exactly as in the scalar case. In the next-to-leading 
approximation, we correct the hard degrees of freedom by their interaction 
with the soft modes. That is, we continue to use the hard thermal loops 
at small momenta, but use at hard momenta the corrections corresponding 



to the diagrams displayed in Fig. 7.5. The resulting approximation to the 



entropy, S = Snla accounts then fully for the orders g and g . But of 
course these expressions are not limited to values of the coupling as small as 
required for the validity of perturbation theory. 

7.6 Some results for QCD 

As an illustration of the quality of the results that are obtained within that 



scheme, we show in Fig. [7^ the entropy of pure SU(3) gauge theory. The 
bands delimiting the various lines in this figure correspond to varying the 
MS renormalization scale p., which defines the renormalized coupling constant 
g{fi), from ft = ttT to AttT. One sees that in contrast to ordinary perturbation 
theory, going from one level of approximation to the next one is indeed a 
small correction. In particular the effects of the soft modes is here a small 
contribution. This is o be contrasted with perturbation theory where the order 
g^ contribution is large for moderate values of the coupling. The comparison 
with the lattice data [^ is quite good down to T > 2.5Tc. 

The quality of the agreement between the self-consistent approximation 
and the lattice data supports the quasiparticle picture of the quark-gluon 
plasma: the dominant effect of the interactions at high temperature seems 
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S/SsB 

1^ 




2.5 



3.5 



4.5 



T/T, 



Fig. 9. The entropy of pure SU(3) gauge theory normalized to the ideal gas entropy 
So- Full lines: Shtl- Dashed-dotted lines: Snla- 2- loop /3- function —» the running 
coupling constant as{jl). The MS renormalization scale: Ji — nT ■ ■ -A-kT. The dark 
grey band: lattice result by Boyd et al (1996) 



to be to change the bare quarks and gluons into massive quasiparticles, with 
small residual interactions between the quasiparticles. It should be empha- 
sized that, in contrast to the approximations based on dimensional reduction, 
the method makes full use of the spectral information on the quasiparticles 
contained in particular in the hard thermal loops. 

The approach is easily extended to finite chemical potential, and the cal- 
culation of the baryonic density can be done using approximations similar to 
those we used for the entropy. Furthermore, from the knowledge of N{ii, T) 
and S{iJL, T) one can reconstruct P(/i, T). Use lattice data to fix the integra- 
tion constant (e.g. P{ji — 0,T)). Such investigations are underway. 
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